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Abstract 



The prospects for the discovery of neutral Higgs bosons (0^ = H^, h^,A^) 
produced with bottom quarks via Higgs decays into muon pairs {pp — > bbcp^ — > 
bbfj,fl + X) at the CERN LHC are investigated in the minimal supersymmet- 
ric model. The complete physics background from the production of bb^ip,, 
bbW~^W~ (including tt) and = g,u,d,s,c in the Standard Model is 

calculated with realistic acceptance cuts. This discovery mode has a simple 
production mechanism from gg — > bbcp^ with its cross section proportional to 
1 / cos^ (3 and could provide an opportunity to measure tan (3 and the bbcf)^ 
couplings. In addition, we compare the associated discovery mode above with 
the inclusive discovery channel pp ^ (fp ^ fifi + X. Promising results are 
found for the CP-odd pseudoscalar (A^) and the heavier CP-even scalar (H^) 
Higgs bosons for tan/? = V2/V1 > 14 and mA,rnH ^ 325 GeV. 
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Typeset using REVT^ 



I. INTRODUCTION 



In the Standard Model (SM), the masses of gauge bosons and fermions are generated by 
one scalar field doublet. After spontaneous symmetry breaking, a neutral CP-even Higgs 
boson {h^sM) remains as a physical particle. 

A supersymmetry (SUSY) between fermions and bosons provides a natural explanation 
of the Higgs mechanism for electroweak symmetry breaking (EWSB) from radiative cor- 
rections in the framework of a grand unified theory and preserves the elementary nature 
of Higgs bosons. The Higgs sector of a supersymmetric theory must contain at least two 
SU (2) doublets |]T| for anomaly cancellation. In the minimal supersymmetric standard model 
(MSSM) 0, the Higgs sector has two doublets (pi and 02 that couple to the ts = —1/2 and 
= -1-1/2 fermions, respectively. After spontaneous symmetry breaking, there remain five 
physical Higgs bosons: a pair of singly charged Higgs bosons if^, two neutral CP-even 
scalars (heavier) and (lighter), and a neutral CP-odd pseudoscalar A^. The Higgs 
potential is constrained by supersymmetry such that all tree-level Higgs boson masses and 
couplings are determined by just two independent parameters, commonly chosen to be the 
mass of the CP-odd pseudoscalar {niA) and the ratio of vacuum expectation values (VEVs) 
of the neutral Higgs fields (tan/? = V2/V1). 

In the Model II of two Higgs doublet models and in the MSSM, the couplings and 
production cross sections of 0° = H^,h^,A^ are enhanced by factors of 1/ cos/5. For 
large tan /3, the rf decay mode [^,0 is a promising discovery channel for the and the 
at the CERN Large Hadron Collider (LHC). It has also been suggested that neutral Higgs 
bosons might be observable via their hb decays [^,0 in a large region of the (m^,tan/3) 
plane, provided that sufficient 6-tagging capability can be achieved. However, simulations 
for the ATLAS detector concluded that detection of the hh channel would be difficult p|. 

The LHC discovery potential of the muon pair channel for neutral Higgs bosons in the 
minimal supersymmetric model was demonstrated by Kao and Stepanov [^|10|! and was 
later confirmed by the ATLAS collaboration 0Jll|. In the minimal supergravity unified 



model [12 , the significance of — >• — > /x/i -|- X is greatly improved at large tan/5 |13 



because and mu become small from the evolution of renormalization group equations 
with large 660° couplings. On the other hand, it is very challenging to search for muon pairs 
from Higgs decays in the Standard Model [0,^. 

Although the /i/i channel has a small branching fraction, this is compensated by an 
improved mass resolution for the muon pairs. For large tan/5, the muon pair discovery 
mode might be the only channel at the LHC that allows precise reconstruction of the A^ 
and the masses. 

In this article, we present the prospects for discovering the MSSM neutral Higgs bosons 
(0° = /i'', A^) produced with bottom quarks via Higgs decays into muon pairs {pp 
hhcj)^ bb^fi + X) at the LHC. We calculate the Higgs signal and complete SM background 
of bbfifi with realistic cuts and compare the discovery potential with that of the inclusive final 
state of pp ^ (j)^ ^ fifL + X. The production cross sections and branching fractions of Higgs 
bosons are discussed in Section II. The dominant physics backgrounds from production of 
bbfifl, and bbW~^W~ are presented in Section HI. The observability of the dimuon signal 
is discussed in Section IV for conservative assumptions about the detector mass resolution. 
Conclusions are drawn in Section V. 
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II. THE PRODUCTION CROSS SECTIONS AND BRANCHING FRACTIONS 



We calculate the cross section at the LHC for pp hhcfP + X (0° = A^) from 



two subprocesses gg — >■ hhcfr ||T6[ and qq bbcj) with the parton distribution functions 
of CTEQ5L ll^. It turns out that the contribution from quark-antiquark annihilation 
{qq — > 660°) is negligible. 

In the MSSM, gluon fusion {gg 0°) is the major source of neutral Higgs bosons 
for tan/3 less than about 4. If tan/? is larger than 7, neutral Higgs bosons in the MSSM 
are dominantly produced from {gg bb(f)^) [T^. Since the Yukawa couplings of (p^bb are 
enhanced by 1/ cos/5, the production rate of neutral Higgs bosons is usually enhanced at 
large tan/?. For rriA larger than about 150 GeV, the couplings of the lighter scalar to 
gauge bosons and fermions become similar to those of the SM Higgs boson. In this case, 
gluon fusion is the major source of the even if tan/? is large. 

In general, the QCD radiative corrections to the subprocess gg 0° are substantial 
||T8|^T1. However, the next to the leading order (NLO) cross section of pp —>■ cfPbb + X via 
gg (fPhb with QCD radiative corrections has been found to be only slightly reduced from 
the leading order (LO) contribution evaluated with the bottom-quark pole mass . 
It is a good approximation to evaluate the NLO cross section of pp (fPbb + X with the 
leading order (LO) contribution multiplied by a K factor[| of 0.8 and using the pole mass 
Mfc = 4.7 GeV. 

The results of Ref. are applied to compute the NLO rates for pp — > 0°. We take 
a K factor of 1.5 for the contribution from gg —>■ H^, hP for tan/3 < 6 and a K factor of 
1.3 for tan/3 > 6. The contribution from gg — ^ is calculated with the LO contribution 
multiplied by a K factor of 1.7 for tan/3 < 6 and a K factor of 1.3 for tan/3 > 6. 

The cross section for pp — > 660° — > bbfifi + X is evaluated with the Higgs production cross 
section a{pp — > 660° + X) multiplied by the branching fraction of the Higgs decay into muon 
pairs 5(0° fifi). With QCD radiative corrections to 0° — ^ 66 ||2^, the branching fraction 
of 0° fifi is about m^/3mf (m^) when the 66 mode dominates Higgs decays, where 3 is a 
color factor of the quarks and mb{'m^) is the running mass at the scale m^. This results in 
the branching fraction for A^ /iji of approximately 3 x 10~^ for = 100 GeV. 



III. THE PHYSICS BACKGROUND 

The dominant physics backgrounds to the final state of bbfifi come from gg bbfifi 
and qq — >• 66/i/i as well as gg — ^ bbW~^W~ and qq bbW^W~ followed by the decays 
of — i> /U^i/^. The Feynman diagrams for gg 66/i/i are shown in Fig. 1. We have 
also considered backgrounds from pp gbfip. + X, pp gbfifi + X, pp ^ gQf^fi + ^, 
pp gqfifi + X, and pp jj/i/i + X, where q = u,d, s,c and j = g,q or q. 



^ The K factor is defined as K = ctnlo/clOj where aio is evaluated with LO parton distribution 
functions and 1-loop evolution of the strong coupling (a^) and ctnlo is evaluated with NLO parton 
distribution functions and 2-loop evolution of a^. This is consistent with the notation of Refs. |19] 
and pi. 
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Our acceptance cuts and efficiencies of 6-tagging and mistagging are the same as those 
of the ATLAS collaboration In each event, two isolated muons are required to have 
PrifJ') > 20 GeV and \ri{^)\ < 2.5. For an integrated luminosity (L) of 30 fb~^, we require 
Prib^i) > 15 GeV and \ri{b,j)\ < 2.5. The 6-tagging efficiency (ej,) is taken to be 60%, 
the probability that a c-jet is mistagged as a 6-jet (ec) is 10% and the probability that any 
other jet is mistagged as a 6-jet (ej) is taken to be 1%. For < 100 GeV, we change the 
requirement to prifJ') > 10 GeV for muons in both the Higgs signal and the background. 

For a higher integrated luminosity of 300 fb~^, we require the same acceptance cuts as 
those for L = 30 fb~^ except Prip,]) > 30 GeV and = 50%. In addition, we require that 
the missing transverse energy ( in each event should be less than 20 GeV for L = 30 fb~^ 
and less than 40 GeV for L = 300 fb~^ to reduce the background from pp bbW^W~ + X 
which receives its major contribution from both real and virtual top quarks pp — > t*i* + X. 

We have employed the programs MADGRAPH p4| and HELAS p5| to evaluate the 
background cross sections of pp — > bbfijl + X,jjfj,p, + X and bbW^W~ + X. In Figure 2, 
we present distributions for the transverse momenta [px) of muons and bottom quarks for 
pp bbA^ —>■ bbfip, + X and for the SM background of pp —>■ bbfip, + X via gg bbfip, 
and qq — ^ bbfifi. It is clear that the muon pt cut is effective in removing most of the SM 
background, while most muons from the Higgs decays pass the pt cut. Most b quarks in 
both the signal and the background do not have large px, therefore, only approximately 10% 
of the bbfifi events survive the requirement of pT{b,j) > 15 GeV for = 100 GeV. 

We have also studied the inclusive discovery channel of — 0° — fifl + X. The 
dominant physics background to the inclusive final state of comes from the Drell-Yan 
process gg — >• Z, 7 — > /i/i . We make the same cuts discussed above and veto events 
with any jet passing the cuts on px and rj to effectively reduce additional backgrounds from 
pp ^ ti ^ bbW^W~ + X and pp jjfJ'fi + X where j = g, u, d, s, c or b. 



IV. THE DISCOVERY POTENTIAL AT THE LHC 



We show the invariant mass distribution of muon pairs in Fig. 3 for pp bbA^ + X ^ 
bbfifL+X with tan (3 = 10 and 50, the SM processes of pp bbfift+X, and contributions from 
pp —>■ bbW^W~ + X with cuts and tagging efficiencies for L = 30 fb~^. Also shown are the 
muon pair invariant mass distributions for the inclusive final state of — > — > fifl + X 
and its dominant background from the Drell-Yan process pp Z,'~f ^ fip, + X. There 
are several interesting aspects to note from this figure: (a) including two b quarks in the 
final state greatly improves the signal to background ratio with realistic efficiencies of b 
tagging and jet mistagging at the ATLAS and the CMS detectors, (b) the SM subprocess 
of gg — s> bbfifi and gg bbfifi make the major contribution to the physics background for 
M^fi < 100 GeV, but gg — > bbW^W~ and gg bbW^W~ become the dominant background 
for higher muon pair invariant mass. 

To study the discovery potential of pp — * 
the background from the SM processes of pp 
where AM„„ = 1.64[(r<^/2.36)2 + a^U/a 



11 



° + X — * bb^fi + X at the LHC, we calculate 
^ bbfifl + X in the mass window of m^±AMf^p, 
is the total width of the Higgs boson, and 
to be 2% of the Higgs boson mass [O]. The 



(Tm is the muon mass resolution. We take cr„ 
CMS mass resolution will be better than 2% of for < 500 GeV [miO|- Therefore, the 
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observability for the muon pair discovery channel at the CMS detector will be better than 
what is shown in this article. 

To study the observability of the inclusive muon discovery mode, we consider the back- 
ground from the Drell-Yan process, gg ^ Z, 7 — fifi, which is the dominant background 



PHTl|. We require (i) prifJ') > 20 GeV and (ii) \ri{fi)\ < 2.5 for both the signal and back- 
ground. 

We define the signal to be observable if the Na lower limit on the signal plus background 



is larger than the corresponding upper limit on the background p6|,^, namely, 



L{as + ab) - NJL{as + (^b) > Lab + NJha. 



which corresponds to 



a., > 



1 + 2JLab/N 



(1) 



(2) 



Here L is the integrated luminosity, o"s is the cross section of the Higgs signal, and ab is the 
background cross section within a bin of width ±AM^^ centered at m^. In this convention, 
N = 2.5 corresponds to a 5a signal. We take the integrated luminosity L to be 30 fb~^ and 
300 fb"^ ini. 

For tan/3 ^ 10, tua and mn are almost degenerate when niA ^ 125 GeV, while and 
rrifi are very close to each other for ^ 125 GeV in the MSSM @|l3. Therefore, when 
computing the discovery reach, we add the cross sections of the and the for < 125 
GeV and those of the A° and the for > 125 GeV. 

The 5a discovery contours for the MSSM Higgs bosons at y/s = 14 TeV for an integrated 
luminosity of L = 30 fb~^ and L = 300 fb~^ are shown in Fig. 4. We have chosen 
^susY = fTT'q = iTT'g = = = 1 TcV. The discovcry region is the part of the parameter 
space above the 5a contour. 

We have presented results for Msusy = 1 TeV. If Msusy is smaller, the discovery region 
of — > /i/i will be slightly reduced for niA <L 250 GeV, because the the Higgs bosons 
can decay into SUSY particles |^ and the branching fraction of 0° — > /x/i is suppressed. For 
niA ^ 125 GeV, the discovery region of /i/i is slightly enlarged for a smaller Msusy, 

but the observable region of /i/i is slightly reduced because the lighter top squarks 

make the and the lighter; also the H^bb coupling is enhanced while the h^bb coupling 
is reduced [1,0]. 



V. CONCLUSIONS 

The muon pair decay mode is a promising channel for the discovery of the neutral Higgs 
bosons in the minimal supersymmetric model at the LHC. The and the should be 
observable in a large region of parameter space with tan (3 > 10. The should be observable 
in a region with niA ~ 125 GeV and tan/? > 5. 

In particular. Fig. 4 shows that the associated final state of 660° bbfift could discover 
the A^ and the at the LHC with an integrated luminosity of 30 fb~^ if ^ 300 GeV. 
At a higher luminosity of 300 fb~^, the discovery region in rriA is not expanded much; the 
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harder pt cut on b quarks reduces the Higgs production cross section, and the higher missing 
transverse energy sUghtly increases the background from bbW^W^ . 

The inclusive final state of 0° njl could also allow the discovery of the and the 
at the LHC with an integrated luminosity of 30 fb~^ if rriA ^ 450 GeV. At a higher 
luminosity of 300 fb~^, the statistical significance of the Higgs signal is improved by a factor 
of \/T0 and the discovery region in tua is significantly extended to ruA ^ 650 GeV for 
tan (3 = 50. 

The discovery of both 0° — > rf and 4>^ — > njl will allow us to understand the Higgs 
Yukawa couphngs with the leptons. The discovery of the associated final state of bbcf)^ 
bb/ip, will provide information about the Yukawa couplings of 660° and an opportunity to 
measure tan f3. Although the discovery region of the fifi mode is smaller than the rf channel, 
the njl channel allows a precise reconstruction of the Higgs boson masses. Furthermore, if 
tuh — TTiA is smaller than the tau mass resolution, it would be difficult to separate them in 
the tau pair discovery mode. The muon discovery channel at the LHC might be employed 
to distinguish the and the A° in multi-Higgs-doublet models for mn — tua shghtly larger 
than the muon mass resolution ((T^) of the CMS and the ATLAS detectors. 
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FIG. 1. Feynman diagrams for gg — > bbmi. 
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Vs"= 14 TeV, Ma = 100 GeV, tan/S = 10 
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FIG. 2. The pT distributions of muons and bottom quarks in ih/GeV for 

pp — > bbA^ bbfijl + X (solid and dot-dashed) with niA = 100 GeV, tan/3 = 10 and 
rriq = rrig = fj, = 1 TeV. We have chosen fii or 6i to be the muon or the b quark with higher 
transverse momentum. Also shown is the contribution from the SM background of pp bb^fl + X 
(dashed and dotted) via gg — > bbufi and qq bbuJl. The pT cuts remove events to the left of the 
vertical dashed line. 
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The invariant mass distribution (da/dM^p) for pp — > bbA^ + X 



bbfifi + X at 



^/s = 14 TeV, as a function of M^^ = niA, for rriq = = /i = 1 TeV and tan/3 = 10 or 50. 
Also shown are the SM backgrounds from pp — > bb^fi + X (dashed), pp — > jj'/U/i + X (dotted), and 
pp bbW~^W^ + X (dot-dashed). The panel on the right shows the invariant mass distribution 
for the inclusive final state pp ^ ^ ^fi + X and the Drell-Yan background (dashed) with 
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MM 



ruA- We have applied the acceptance cuts as well as the efficiencies of b tagging and 



mistagging for L = 30 fb 
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MSSM, MsusY = 1 TeV 

(a) L = 30 fb"^ (b) L = 300 fb"^ 
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FIG. 4. The 5a contours at the LHC for an integrated luminosity (L) of 30 fb^^ and 300 fb~^ 
in the niA versus tan /? plane. The signal includes (j)^ = and for niA < 125 GeV, but (p^ = 
and for > 125 GeV. The discovery region is the part of the parameter space above the 5a 
contour. 
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